Interference of two electrons entering a superconductor 

F. W. J. Hekking and Yu.V. Nazarov 
Institut fur Theoretische Festkorperphysik, Universitat Karlsruhe, Postfach 6980, 7500 Karlsruhe, 
FRG 



The subgap conductivity of a normal-superconductor (NS) tunnel junction is 
thought to be due to tunneling of two electrons. There is a strong interference be- 
tween these two electrons, originating from the spatial phase coherence in the normal 
metal at a mesoscopic length scale and the intrinsic coherence of the superconductor. 
We evaluated the interference effect on the transport through an NS junction. We 
propose the layouts to observe drastic Aharonov-Bohm and Josephson effects. 
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Quantum phase coherence in solids manifests itself basically in two ways. First, there is 
an intrinsic coherence in the superconducting state. In superconductors, the phase is indeed 
a macroscopic variable. This can be observed in a variety of interference experiments, for 
example, with tunnel junctions Second, even in a normal metal electrons are coherent 
at a mesoscopic length scale. Interference between the electrons in a normal metal gives rise 
to a set of phenomena which constitutes the subject of a new branch of condensed matter 
physics i|. 

At the NS interface between a normal metal and a superconductor, these two sources of 
coherence may interplay. Very recent experiments j|] show how interesting such an interplay 
may be. These experiments were performed with NS boundaries of a high transparency. The 
phenomenon of Andreev reflection |5j seems to be responsible for the peculiarities observed. 

In the present paper we focus on the opposite case of a tunnel NS interface. It is well- 
known that an electron with energy less than A, A being the superconducting energy gap, 
can not tunnel to the superconductor and therefore the transport through the junction is 
strongly suppressed at voltages below the gap [|TJ. On the other hand, two electrons can 
enter the superconductor converting into a Cooper pair since this process costs no energy. 
Such two-electron tunneling gfj] determines the subgap conductivity of the junction. It 
has been discussed f| that the rate of this two-electron process is often determined by the 
interference of the electron waves on a space scale given by the coherence length, either in 
the normal or the superconducting metal. That motivates us to explore how agents which 
act on the phase will influence the subgap conductivity. Indeed we find the conditions under 
which a pronounced Aharonov-Bohm and Josephson effect can be observed. 

Since we consider the low- voltage subgap conductivity, we assume that T, eV <C A. 
Under these condictions the coherence length = ^JhD/ksT, where D is the diffusion 
constant, is much larger than the one in the superconductor, ■JhD/A. Therefore we con- 
centrate on the interference in the normal metal. 

The NS interface is described by the Hamiltonian H = Hn + Hs + Ht, where Hn 
and Hs refer to the normal and the superconducting electrode, respectively. The tunnel 
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Hamiltonian Ht is given by the usual form Ht = Ek,p,u t^p^L J*?^ + t^ p b pa dk^. Here, 
operators a^, b p>a correspond to the normal and the superconducting electrode, ik P are 
the tunnel matrix elements which we take to be spin-independent; the sum is taken over 
momenta k, p and spin a =1,1- Second order perturbation theory in Ht yields the lowest 
order contribution to the amplitude of two-electron tunneling. The a operators appearing 
in this amplitude remove two electrons from the normal metal electrode with energy and 
£/./. The amplitude thus consists of a sum over intermediate states in the superconductor 
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The matrix elements between \N) and |7V — 2) connect states differing by two electrons. In 
coordinate representation they can be expressed in terms of the usual anomalous Green's 
function F^(r[, r' 2 ; u). Since we assume that T, eV A, we consider only (,k,£,k' *C A and 
find 
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f(n,r' 1 )f(r 2 y 2 )Fl(r' i y 2 ;u, = 0) (2) 

where ipk{ r ) denotes an eigenfunction of an electron in the normal metal; primed space 
arguments refer to the superconductor. In a disordered material, ipk( r ) and Pl,(r[, r' 2 \ uS) are 
in general complicated functions depending on the realisation of the disorder. Since we are 
interested only in the interference occurring in the normal metal, we perform an average of 
(0) over states in the superconductor. This may be done along the lines of flTI ; the product 



of two tunnel amplitudes will give the normal state conductance g(r) of the tunnel interface 
per unit area such that the total conductance Gt = J d 2 rg(r). The result reads 



A W = ~ 2 — I d 2 rg(r)Mr)Mr)e^ r) (3) 
e z u N J 



where 0(r) is the phase of the superconducting condensate and vn the density of states in 
the normal metal. The rate for two-electron tunneling as a function of the applied bias 



voltage V is obtained by applying Fermi's Golden Rule. To obtain the current we have to 
sum the tunnel rates in both directions. As a result we find 

I(V) = — — / d 2 r 1 d 2 r 2 g(r 1 )g(r 2 ) expi^r^ -0(r 2 ))x 

J du{f{u/2 - eV) - f(u/2 + eV)}{pC(r u r 2 ) + P? u (n, r 2 )} (4) 

Here / is the Fermi distribution for electrons in the normal metal. Eq. @ is the central result 
of our paper, which clearly shows the interplay between phase coherence in a superconductor 
and a normal metal. The intrinsic coherence of the superconductor is reflected by the 
appearance of the phase difference <f>(rx) — <f){r 2 ). In the normal metal, the two incoming 
electrons undergo many elastic scattering events in the junction region before they tunnel 
through the NS interface, leading to interference on a length scale given by Lt ||. These 
interference effects have been taken into account by averaging the rate in the standard way || 
over possible scattering events. The result (f|) therefore contains the sum of two Cooperon 
contributions Pj(r, r'), which obey the equation [11 



{-%D{V - i2nA(r)/%) 2 - ico}P^'(r, r') = 5{r - r') (5) 

where A is the vector potential and $o — hc/2e the flux quantum. From this equation it 
is clear that the result does not only depend on properties of the junction (via Gt), but 
also on its surroundings over a distance Lt, due to the interference occurring on this length 
scale. 

As a simple example we calculate first the subgap conductance corresponding to a layout 
where a semi-infinite normal wire of thickness d <C Lt is connected to a superconducting 
electrode by a tunnel junction. In this case, the only contribution from the spatial integra- 
tions in (^) originates from the tunnel junction at r\ = r 2 = 0. The solution of (Bl) for a 
wire yields (0, 0) = l/dy/-iuhD. It leads to 

G wire = 8vr 5 / 2 (l - 2V2)C{-l/2)R cor G 2 T w 53.8^^ (6) 

where R cor = LT/{e 2 u^Dd) is the resistance of the wire per correlation length Lt- 



Interference effects in a mesoscopic system threaded by a magnetic flux lead to the 
Aharonov-Bohm effect: the total resistance depends periodically on the applied flux ||12|| . A 
similar effect can be observed with the layout depicted in Fig. la, where a small loop with 
circumference L is inserted into a wire at distance I of the junction. The resistance of the 
loop is denoted by Rl, Ri is the resistance of the piece of wire between loop and junction. 
The loop is threaded by a magnetic flux $. The conductance (f|) at zero temperature 
{Rl,Ri < Rcor) reads 

G loop = 4tt 2 G 2 t + . 2 ^.. | (7) 
[ snr7r$/$oJ 

This result is plotted in Fig. 2 (upper curve). It shows the usual h/2e periodicity |12| related 
to the Cooperon. Moreover it diverges each time when $ = $o- At finite temperature 
the result (|6]) restricts the maximal conductance by the value 53.8R cor G^. But even at 
zero temperature, when R cor diverges, the penetration of magnetic flux in the wires that 
constitute the loop. This penetration leads to a shift in the Cooperon energy uo — > to + 
ia 2 (hD/4:L 2 )(2n^>/^ ) 2 where a = S W i re /Si oop is the ratio of the area of the wire and the 
loop. As a result we obtain instead of (0): 

_ 8ir 2 G 2 T R L (l + 2(l/L) sm 2 7r1>/%) 2 

loop (7ra<l>/2$o)(cos27r$/$o + x ) + (2 sin 2 tt$/$ )(1 + 2(1/ L) sin 2 vr$/$ ) 1 j 

Thus the divergencies are removed as can be seen in Fig. 2 (lower curve), where the result 
is plotted, taking l/L = 0.5, and a = 0.1. We note that the flux-dependence presented here 
will be absent when the superconductor is in the normal state. 

We now turn to a different geometry (Fig. lb), where instead of a ring a fork is attached 
to the wire, such that we have two tunnel junctions to the superconductor at different 
positions r\ and T2- The subgap conductance will be determined not only by the flux 
threading the closed area between the fork and the superconductor, but also by the magnetic 
field distribution in the superconductor. Let us consider the curves Cn and Cs connecting 
the junctions 1 and 2 in the normal metal and the superconductor, respectively. The effect 
will be governed by the phase 9 = 0(r x ) — 0(r 2 ) + (2e/hc) J Cn A.dx. In order to obtain a 
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gauge-invariant expression for 9 we use the relation V0 — (2e/hc)A = p s /h, where p s is the 
momentum of the superconducting condensate. In this way we arrive at 9 — 27r$/$ + 
fc s (Ps/h)-dx, with $ the flux penetrating the closed loop formed by Cn + C$- This result 
does not depend on the choice of Cs- It reflects the dependence of 9 on the penetration of 
the magnetic field as well as on the vortex positions in the superconductor. The effect can 
be used to monitor these positions. Introducing the conductances Gi,G?2 of the junctions 
we obtain 

(9) 

when Lt is larger than the size of the fork. The conductance Gf or k thus combines the phase 
coherence in the normal and the superconducting metal. 

If we consider tunneling to a superconductor of finite size, the transport will be influenced 
by charging effects, if the total capacitance C of the superconducting island is small enough, 
such that the charging energy E c = e 2 /2C is of the order of A 0,|lj. We will restrict 
ourselves to the case A > E c . In this case the transport to the superconductor will be 
due to two-electron processes. We will present a simple relation between the rate for these 
processes and the current which would flow in the absence of charging effects. The transfer of 
two electrons to the superconducting island increases its electrostatic energy by an amount 
E2ei- This energy can be changed with the help of an additional potential V g applied to 
the island; we assume that E 2e i <C E c , A. The first electron entering the superconductor 
as a quasiparticle then increases the electrostatic energy of the island by an amount E c . 
The energy of the intermediate states in the superconductor will therefore be shifted by this 
amount, leading to a different amplitude v4 k | k /j = (F(E C )/ F(0))Ay.^>^, with 

F(E) = ; 4A arctan J - + f (10) 
V ; VA 2 - E 2 \A-E v ' 

The rate is then given by 2ef = (F{E C )/ F{0)) 2 I(E 2e i/2e). 

Finally we investigate a Josephson-like effect, which occurs in the geometry depicted 
in Fig. lc. In this layout a normal fork is connected to two different superconductors, 
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G 2 + G 2 . + 2G 1 G 2 cosfl 



to which a small voltage difference eVs is applied. The phase difference between the two 
superconductors and hence between the junctions at the extensions of the fork will increase 
linearly with time: 9 = eVst/h. Substituting this phase difference into Eq. (0) we obtain 



G T = 53.8/L 



G\ + G\ + 2G X G 2 cos(eV s t/h) 



(11) 



Thus the conductance oscillates with a frequency ujj = eVs/h. The modulation is of the 
order of the conductance itself. 



In conclusion we studied the effect of interference on the subgap conductivity of an NS 
tunnel junction. Transport is determined by the transfer of electrons in pairs from the normal 
metal to the superconductor. At low temperatures, interference between the two electrons 
occurs in the normal metal over a longer length scale than in the superconductor. Therefore, 
the subgap conductance is determined not only by properties of the tunnel interface, but 
also on the layout on the normal side near the interface over a distance Lj-. These novel 
interference effects can be made visible by influencing the electron phase, e.g. with the help 
of the Aharonov-Bohm effect or the Josephson effect. We discuss these effects for various 
layouts of practical interest, and present results for the sugap conductance. 
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FIGURES 

FIG. 1. Three geometries discussed in the text: (a) 1-dimensional normal wire containing a 
loop connected to a superconducting electrode by a single junction; (b) fork geometry connected to 
the superconductor by two junctions; (c) fork geometry connected to two different superconductors. 

FIG. 2. Subgap conductance at zero temperature for the geometry of Fig. la, as a function of 
flux. Curves correspond to Eqs. (0) (upper curve) and (|8|) (lower curve) in the text. 
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